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Abstract 

The Casimir mutual free energy F for a system of two dielectric 
concentric nonmagnetic spherical bodies is calculated, at arbitrary 
temperatures. The present paper is a continuation of an earlier inves- 
tigation [Phys. Rev. E 63, 051101 (2001)], in which F was evaluated 
in full only for the case of ideal metals (refractive index n = oo). Here, 
analogous results are presented for dielectrics, for some chosen values 
of n. Our basic calculational method stems from quantum statisti- 
cal mechanics. The Debye expansions for the Riccati-Bessel functions 
when carried out to a high order are found to be very useful in practice 
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(thereby overflow/underflow problems are easily avoided), and also to 
give accurate results even for the lowest values of I down to Z = 1. 
Another virtue of the Debye expansions is that the limiting case of 
metals becomes quite amenable to an analytical treatment in spheri- 
cal geometry. We first discuss the zero-frequency TE mode problem 
from a mathematical viewpoint and then, as a physical input, invoke 
the actual dispersion relations. The result of our analysis, based upon 
the adoption of the Drudc dispersion relation at low frequencies, is 
that the zero-frequency TE mode does not contribute for a real metal. 
Accordingly, F turns out in this case to be only one half of the con- 
ventional value at high temperatures. The applicability of the Drude 
model in this context has however been questioned recently, and we 
do not aim at a complete discussion of this issue here. Existing exper- 
iments are low-temperature experiments, and are so far not accurate 
enough to distinguish between the different predictions. We also cal- 
culate explicitly the contribution from the zero-frequency mode for a 
dielectric. For a dielectric, this zero-frequency problem is absent. 
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1 Introduction 



In the Casimir world, it is desirable to consider geometrical configurations 
that are amenable to an analytical treatment and at the same time nontrivial 
enough to elucidate the physically important properties. The following con- 
figuration turns out to satisfy these two criteria (cf. Fig. 1): there are two 
spherical bodies present with concentric surfaces at r = a and r = b, with a 
vacuum region in between. We shall consider the free energy F{T) due to the 
mutual interaction between the two bodies. We gave earlier an analysis of 
this problem 0], with the use of quantum statistical methods as well as field 
theoretical methods. Whereas the general formalism in was valid for ar- 
bitrary (equal) permittivities e in the two dielectric regions r < a and r > b, 
the explicit evaluation of F{T) for various values of T and widths d = (b — a) 
was made for the case of perfectly conducting walls only, corresponding to 
e ^ oo. Our purpose with the present paper is to extend these previous 
considerations to cover the case of general values of the permittivity. As to 
our knowledge such a calculation has not been undertaken before, although 
there are similarities with the theory given by Kleinert some years ago [Q]. 
We will henceforth assume, as in that the two media are nonmagnetic. 
A brief account of the essentials of the present theory was recently given in 
Ref. [0]. We will have to repeat some of the formalism below because of 
readability. 

One lesson from the calculation in |l|] was that the power of the quan- 
tum statistical method is remarkably strong when applied to the rather 
demanding case of general e. The most central formula in our context is 
the statistically-derived Eq. (40) in it gives the value of PF = F/T 
for arbitrary values of temperature, width, and e. Whereas this equation 
was given in terms of a very compact notation in it will be convenient 
here to rewrite it slightly. Let m G {—oo, oo) be an integer correspond- 
ing to Matsubara frequencies K = 27im/(3; let n = ^/e be the refractive 
index of the two media lying at r < a and r > b, and let si{x), ei{x) 
be Riccati-Bessel functions with imaginary argument defined according to 
si{x) = (ttx/2Y^'^I^{x), ei{x) = {2x/7fY^'^K^{x), so that their Wronskian 
becomes W{si,ei} = —1. Here u = I + 1/2, and Jj,, K^, are modified Bessel 
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functions. We write the formula as 

oo ' oo 

/^^ = E E(2^ + - A^) + ln(l - XD], (1) 

m=0 1=1 

where the prime on the summation sign means that the m = term is taken 
with half weight. The two eigenvalues Af*^ and Xf^ in Eq. (1) correspond to 
the transverse magnetic and the transverse electric modes. (In the notation 
of Ref. Ill, Xei = A™, A; = Af^; we find it here useful to emphasize 
the physical nature of the two modes.) For later use we will write these 
eigenvalues as ratios. First, 

\TA4 _ fih /o^ 
J3J4 

where 

/i = ns'i{x)si{nx) — Si{x)s'i{nx), 
f2 = ne'i{y)ei{ny) - ei{y)e'i{ny), 
/g = ne'i{x)si{nx) - ei{x)s'i{nx), 

h = nei{ny)s'i{y) - e[{ny)si{y), (3) 
x and y being the nondimensional frequencies 

X = 2nma/P, y = 2Txnih/(3. (4) 

We put h = c = kB = 1] /? = 1/T is the inverse temperature. It should be 
noted that, in contradistinction to the formalism in the primes in Eqs. (3) 
mean derivatives with respect to the whole argument. 
Next, the TE eigenvalues are written as 

xr = ^, (5) 

where 

gi = s'i{x)si{nx) — nsi{x)s'i{nx), 
92 = e[{y)ei{ny) - nei{y)e\{ny), 
= e'i{x)si{nx) — nei{x)s'i{nx), 
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Qa = ei{ny)s'i{y) - ne[{ny)si{y). 



(6) 



The case of metaUic walls, n = ^/e — > oo, leads to a delicate two-limit 
problem as regards the contribution from zero Matsubara frequency, m = 0. 
The conventional way to proceed when handling this problem within the 
framework of nondispersive theory, has been to take the limits in the following 
order: 

(i) Set first e = oo; 

(ii) take then the limit m —>■ 0. 

This way of taking the limits was advocated already in the 1978 paper 
of Schwinger, DeRaad, and Milton , and the same procedure was followed 
in Sec. VII of our previous paper |T|]. The method implies inserting the 
small-argument approximations for the Riccati-Bessel functions into Eq. (1), 
resulting in the following free energy expression: 



m=0 1=1 



sijx) ei{y) 
ei{,x) si{y) 



1 



e'lix) s'liy) 



(7) 



which is in agreement with Eq. (68) in If we next let a; — > 0, y — > 
observing the same low-argument approximations, we obtain as contribution 
from 171 = 0: 



ppconv^^ ^ oo^ ^ = 0) = ^(2/ + 1) In 



1=1 



2i+i 



(8) 



again in agreement with [jT|, Eq. (79). This is the conventional result. Both 
the two electromagnetic modes are in this way found to contribute equally 
to the sum in Eq. (8). 

However, a discussion has recently arisen as to whether this recipe for 
dealing with the m = term in the TE mode is really correct. The problem 
becomes most acute in the high T regime, but is present at moderate and low 
temperatures also. We may refer to the paper of Bostrom and Sernelius 
questioning this point, and the subsequent comment of Lamoreaux [0. What 
has been most welcome in recent years are the accurate experiments on the 
Casimir force, due to Lamoreaux [0], Mohideen et al. P, |10|, [ill, and Bressi 



et a/.[r^. By means of these experiments it becomes much easier to formulate 
a sound theory. Several theoretical papers have appeared, discussing various 
facets of the experiments [1^, 0, ^7^. An extensive recent review 
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has been given by Bordag et al. [|T^]. We also mention several other related 
papers |19, 20, 21 1, of a more general nature, though also these are being 
concerned with finite temperature effects in a Casimir context. 

One of the purposes of the present work is to analyze how the m = case 
works out for case of the spherical geometry of Fig. 1. It turns out that the 
formalism becomes actually quite manageable. Use of the Debye expansion 
for the Riccati-Bessel functions is an essential element in our analysis, and 
it implies as an additional bonus that the overflow and underflow problems 
that so often plague calculations of this sort, are easily abandoned. 

Our analysis obviously requires machine evaluation, as we will carry out 
the Debye expansion to the 18th order in the quantity 9 deflned in Sec. 
2.1 below. Then the numerical accuracy becomes quite satisfactory for all 
practical purposes, even for the lowest values of / down to / = 1. Up to 
about one million terms in the series will be summed. The limiting case 
of metals will be handled in a physical rather than a mathematical way by 
adopting the physically preferable Drude dispersive model as input at low 
frequencies. On basis of the Drude dispersion relation, we are quite naturally 
led to the conclusion that there is no contribution from the m = TE mode 
to the free energy F in the case of inflnite conductivity. This implies that 
the conventional expression for F for metals has to be multiplied by one 
half. (This conclusion is actually in agreement with the outcome of the 
statistical-mechanical considerations in Sec. Ill in [|[].) We show graphically 
several results for how F varies with temperature and width, both for the 
ideal metallic case (in which the zero-frequency mode is counted twice) and 
for the dielectric case. Finally, we calculate the magnitude of the m = 
contribution to F for a dielectric, as a function of temperature, and compare 
the result with the total value of F. Also, the mutual internal energy E itself 
is briefly discussed. 

We do not in the present paper aim at resolving the issue about the m = 
term for the TE mode for real metals. We make however some estimates in 
the discussion in Sec. 5, item 6, choosing aluminum as a concrete example. 
As for experiments, the atomic force microscope experiment of Mohideen 



and Roy ||^, and that of Harris et al. [jTO[, achieved an accurary of about 
one per cent. These are essentially low-temperature experiments, where the 
influence from the m = TE term is small (at T = the m = is completely 
negligible since the sum over discrete Matsubara frequencies is replaced by 
an integral over imaginary frequencies). The all-over temperature correction 
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at room temperature is predicted to be of the same one per cent accuracy 
. The single m = temperature term which we discuss, is not singled 
out experimentally under these circumstances. 



2 Numerical considerations 

We now define the nondimensional temperature: 



implying that x = mt. It turns out numerically that the conventional 
uniform asymptotic expansions of the Riccati-Bessel functions, which are 
used often to low orders when dealing with rough approximations, become 
quite accurate if the polynomial parts of the expansions are expanded to 
high order. This makes the evaluation in the present case straightforward 
in principle: the polynomial parts, which generally turn out to be about 
unity in magnitude, can easily be handled on a computer. The remaining 
parts of the Bessel functions, which are simple exponentials, can be dealt 
with analytically. Actually, what are needed in practical calculations, are 
fractions between Bessel functions. This is the way in which the overfiow or 
underflow problems are avoided. Overflow/underflow problems would easily 
arise if we instead chose to take the whole Bessel function directly from the 
computer library. The remaining numerical evaluation is not quite trivial, 
though; especially in the case of a narrow slit the number of necessary terms 
turns out to be quite large, about 10^, as we mentioned above. 

We start by presenting our expanded version of the Debye formalism. 

2.1 The Debye expansions 

Let us write the Debye expansions in the form 




(10) 




(11) 
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Si[X) 



ei[x) 



^Jl±^Ml!le''^^^^C[e{x)], 
e-"^^-^ D[e{x)]. 



2 ^z{x) 

[1 + ^2 (3,)] 1/4 



(12) 



z{x) 



Here u = l + 1/2, / = 1,2, 

r]{x) = 



Z[X] 



e{x] 



+ In 



x/v, 9{x) 
z{x) 



1 + i/e{x) 



(13) 

[l + z^{x)]-'/\ and 
(14) 



is the same as the symbol t in Ref. [^). There occur four polynomials 
A{e),B{e),C{e),D{e), which are found to be of order unity. In Ref. |2 



the first two of them, A{6) and B{6), are expanded to order 6^"^, whereas 
C{9) and D{9) are expanded to order 9^. In Ref. |2^ we expanded all the 
polynomials to order 9^^. These expansions, which will not be reproduced 
here, are found to be easily handled on a computer. The polynomials possess 
the important property that they go to unity when 9 goes to zero. 

The factors in Eqs. (10)-(13) that can take extreme values, are the expo- 
nentials. They are easily dealt with analytically. 

It is now convenient to calculate the following ratios between the functions 
defined in Eqs. (3): 



h ^ _1 2.,(.) n'lC[9{x)] - A[9{x)]C[9{nx)]/A[9{nx)] 
h 2 n^-iD[9{x)]+B[9{x)\C[9{nx)\/A[9{nx)y 

.2.,iy) n^D[9iy)]- B[9iy)]D[9iny)]/B[9iny)] 
h n^C[9{y)] + A[9{y)]D[9{ny)]/ B[9{ny)] ' 

where 7 and 5 are the coefficients 
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1 +Z2(x) 

^ 1 + z'^{nx) 
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^ 1 + z'^iny) 



(15) 
(16) 

(17) 



Similarly 



31 = _ 1 -2..(.) lC[9ix)] - A[9ix)]C[9inx)]/A[9inx)] 
gs 2 'yD[9{x)]+B[9{x)]C[9{nx)]/A[9{nx)\ 



(18) 
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91 = _n,-2u,iy) SD[9iy)]- B[eiy)]D[e{ny)]/B[e{ny)] 

g, 6C[e{y)]+A[e{y)]D[e{ny)]/B[e{ny)Y ^ ^ 

Now the eigenvalues Xf^ and Xf^ are calculable from Eqs. (2) and (5), 
with use of the expressions (15)-(19) in which the ^-expansions for the four 
polynomials are taken from Ref. As in our previous paper we made 
use of standard FORTRAN routines throughout. 



2.2 Calculated results for dielectrics 

On a logarithmic plot with base 10, Fig. 2 shows how logiQ{—f3Ft) varies 
with d/a for various values of t when the medium is dilute, n = 1.1. The 
figure is to be compared with the corresponding Fig. 1 in [|l|]. As expected, 
the magnitude of the mutual free energy is much less for a dilute medium 
than it is in the case of ideal metallic walls (n = oo). For instance, when 
d/a = 0.2, t = 1, for n = 1.1 we see that \F\ has only about 0.1 % of the 
value it has for an ideal metal. The various curves in Fig. 2 tend to overlap 
at low temperatures. Thus the curve calculated for t = turns out to be 
indistinguishable from the curve calculated for low temperatures up to t = 1. 
The curves in Fig. 2 are most useful for the case of low temperatures. 

Figure 3 shows how logio(— varies with d/a. This representation is 
convenient for the case of high temperatures, since the curves for high t tend 
to overlap. 

Figure 4 shows the representation in the form that is probably the most 
instructive one, namely log^Q^—PFt) as a function of log^Q/f:. It shows clearly 
how there is a low-temperature plateau, extending up to a region lying some- 
where between 1 and 2 in the cases shown. For higher values of t, there is a 
gradual change into the region where F varies linearly with t. 

We calculated analogous figures for other values of n also, with results as 
one would expect: the influence of the medium becomes strengthened when 
n becomes greater. Figure 5 shows, as an example, the analogue of Fig. 4 in 
the case of n = 2. For instance, when t = 1, d/a = 0.01, the magnitude \F\ 
is about 50 times as large when n = 2 as when n = 1.1. 

Figure 6 shows the analogous variation of F for the case of an ideal metal 
(i.e., n = oo for all uj, the m = mode counted twice). This figure is 
reproduced from Fig. 3 in [|l|; it is included here both for the purpose of 
comparison, and also to correct the labeling on two of the curves in our 
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previous Fig. 3. We remind ourselves that the order of taking the hmits in 
[||] was in accordance with the prescriptions (i) and (ii) mentioned earher, 
above Eq. (7) in Sec. 1 (Refs. 0, |T[). 

Generally, we found the asymptotic Debye expansions to be useful for 
X > 10 and/or I > 9. Then, an accuracy of 8 digits for the individual terms 
was achieved. Below these limits for x and I, we employed the machine- 
generated Bessel functions. For small values of d/a and t, slow convergence 
was observed. The summation of the series thus became rather demanding. 
For instance, when d/a = 0.05, t = 0.01 about 1.1 million terms were needed, 
if we truncated the summation at e = 10~^ (here e means the ratio between 
a general term in the series and the sum). The sum itself is however accurate 
only up to 4 or 5 digits. 

An important result was that even for low values of /, the asymptotic 
series gave very good results. One reason for this is the high-order expansions 
used for the polynomials A, B,C, D. Most probably, the Debye expansions 
(at least when carried out to order 6^^) can be used for all x and /, for all 
practical purposes. 

3 The limiting case of a metal 
3.1 The nondispersive case 

Although it would seem most natural to discuss the case of a metal on the 
basis of a parallel-plates configuration, let us analyze here how the formalism 
behaves in the ideal-conductor limit when the spherical geometry of Fig. 1 is 
given. It will actually turn out that the Debye expansion is very useful also 
in this case. As we treated this topic in reasonable detail in Ref. [§, we need 
only to be brief here. 

We assume first that the medium is nondispersive. The formal limit that 
we have to take, is thus e oo. Let us categorize how to take the two actual 
limits: we let option A mean first taking the refractive index n = ^/e — > oo, 
thereafter taking the Matsubara frequency m — >■ 0. Option B reverses the 
succession of the limits on n and m. 

Consider first the TM mode, employing option A. When n (X), 6{x) is 
finite, while 6{nx) 0. Thus all polynomials {A, B, C, D}[6{x)] at argument 
6{x) are finite, while {A, B,C, D}[6{nx)] 1. Observing that and rP5 
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are proportional to n for large n we get, when taking the limit m ^ 0, the 
following expression for the m = contribution to the TM free energy: 




1=1 



a 



b 



) 



2Z+1 



(20) 



Following instead option B we find precisely the same expression as in 
Eq. (20). The m = TM free energy is thus robust with respect to the 
choice between options A and B. This is actually what we would expect 
physically: the TM mode means that the magnetic field is transverse to the 
radius vector r at r — a,b; this is the natural electromagnetic boundary 
condition at perfectly conducting surfaces. 

Consider then the TE mode, employing option A. The difference from 
the preceding case lies in the sensitivity of \J^{n oo) with respect to the 
coefficients 7 and 5. Prom Eq. (17) we get 7 — >^ 0, 5 — > implying that, in 
the hmit m ^ 0, Xj^ (0/6)^'+^ It follows that the TE contribution to 
the m = free energy is the same as given by Eq. (20). 

Employing instead option B we obtain 7 ^ 1, (5 ^ 1, resulting in Af^ 
when m — > 0. Consequently 



Option B gives accordingly one half of the conventional result of Eq. (8) for 
the total free energy in the ideal conductor limit. 

The immediate question is now: which of the two options is correct? 
We cannot decide upon this only by investigating how the mathematical 
formalism behaves in the limiting cases; we have to bring physics into the 
consideration. That means, we have to consider a physically appropriate 
dispersion relation in the limit of low frequencies. That is the topic of the 
next subsection. 

3.2 The dispersive case 

Let now Cj denote the frequency along the imaginary frequency axis.There are 
mainly two actual dispersion relations on the market. The first corresponds 
to the plasma model of the dielectric: 



B : (3F^^{m = 0) = 0. 



(21) 



e{iLj) = 1 + - 



9 ' 

1-^ 



(22) 
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ujp being the plasma frequency. As mentioned by Landau and Lifshitz (Sec. 
78 in [^]), the range of frequencies over which this formula is applicable 
begins, in practice, at the far ultraviolet for light elements and at the X-ray 
region for heavier elements. Let us for convenience rewrite the coefficients 
(17) as 



7 



1 + {uja/uf 



6 



1 + (cub/u)^ 
\l + {ntbh/uY 



(23) 



When uj — 

n{iCj)u}a/ v 
instance, lo, 



Up, which means that 



0, it follows from Eq. (22) that n{iCj)u) 

Xpjv where, in dimensional units, Xp = Upa/c. Taking, for 
~ 3 X 10^^ s~^ and a ~ 1 cm we get Xp ~ 10^. In practice, the 
most significant values of I are much lower than this. We can thus assume 
that Xp/u ^ 1 in Eq. (23), so that in practice 7 — > 0, 6-^0. That is, we 
recover in this way option A. In conclusion, the use of the plasma dispersion 
relation (38), to a good approximation, leads to the conventional result (8) 
for the m = total free energy for a metal. 

Consider next the Drude model for the dielectric, corresponding to 



:{iuj) = 1 + 



uj{(2i + 7) ' 



(24) 



7 being the relaxation frequency. According to this relation n(ia))a) 
when Co 0, implying that 7-^1, 6^1 according to Eq. (23). That is, we 
recover option B. The total m = free energy for a metal is thus according 
to the Drude model predicted to be one half of the expression (12). 

When deciding between these dispersion relations, we expect that relation 
(24) is physically correct in the limit when uj ^ 0. On general grounds 
the permittivity has to be inversely proportional to the frequency at low 
frequencies; cf. Sec. 77 in Explicitly, e^u) ia/u, or e^iu) = cr/u, 
where a is the conductivity. This is a result following directly from Maxwell's 
equations. The Drude model satifies this requirement. Thus both the Drude 
model (and, as we have seen, statistical mechanical methods), support the 
option B above. The plasma model, Eq.(22), as we have noted, is appropriate 
only for the higher frequencies. 
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4 Calculation of the m=0 contribution to F 
for a dielectric 



The delicate m — problem in the limiting case of a metal accentuates the 
following question: how large is the m — contribution to the free energy 
in the general case, for a dielectric? As the last point in our paper we shall 
calculate this effect, for a given value of n, and show the result graphically in 
a typical example. This point appears to be of physical interest, and with the 
above formalism the calculation can be easily effectuated. We now return 
to nondispersive theory again, so that n is taken to be a constant. Since 
the m — case does not contribute to the TE mode at all, for any finite 
value of n, our present discussion has no bearing on the topic discussed in 
the previous section. 

Let us first summarize, from a physical point of view, how the free energy 
is distributed over the various frequencies for various values of the tempera- 
ture. When T — > the Matsubara frequencies are closely spaced (at T — 
the Matsubara summation being replaceable by an integral), and a large num- 
ber of eigenfrequencies contributes to F. The contribution from the lowest 
term m = is insignificant. When T increases, the number of contributing 
Matsubara terms gradually becomes smaller and the frequencies gathered at 
the lower end of the spectrum until finally, at T ^ oo, the term m — 
dominates completely (this is the classical limit). How this gradual change 
actually occurs, as a function of T, for a given relative slit width d/a, can 
however only by found by an explicit calculation. 

We recall that for a given geometry there are still three quantities to be 
contemplated, namely {n,m,t}. Let us fix the value of n, and look for the 
contribution to F from m = 0, as a function of t. From Eq. (1) we have, for 
an arbitrary temperature. 



(as noted, Xf^ does not contribute for a dielectric). We define Y as the ratio 
between F{m — 0) and the expression (1) for the full free energy: 



pF{m 



0) = ^E(2/+l)ln(l-An 



(25) 



1=1 



F{m = 0) 
F 



(26) 
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For given (i/a, Y thus becomes a function only of t. Figure 7 shows Y versus 
t for various values of d/a, for the case when n = 1.1. The curves behave 
qualitatively as we would expect; the contribution from m = goes to zero at 
very low temperatures, and goes towards unity at high t. Analogous curves 
for other values of n behave similarly; thus the curves calculated for n = 2.0 
turn out to be practically indistinguishable from those in Fig. 7. 

One additional conclusion to be drawn from Fig. 7 is that the less the 
value of d/a, the less becomes the importance of the m = term. It is worth 
noticing that this is a result that can be understood physically: when the slit 
is narrow, as assumed in Fig. 7, we can approximately regard the system as 
a conventional two-plates system. For the latter geometry, it is known that 
the classicality condition can be written as dT ^ 1, where d is the distance 
between the plates (cf. Sec. 82 in [^). When d decreases the system 
thus becomes more and more a quantum-mechanical system, necessitating 
an increasing large region of frequencies determining the value of F. The 
relative importance of the low frequencies, in particular that of m = 0, thus 
has to diminish, in accordance with the figure. 

5 Conclusive remarks 

Let us summarize our work, and supply the above with some further remarks. 

1. The Debye expansion procedure is almost surprisingly effective. When 
carried out to sufficiently high order in the parameter 6 - order 18 in the 
present paper - the accuracy becomes fully satisfactory for all practical pur- 
poses for all values of /, even down to the lowest value / = 1. Moreover 
the formalism becomes straightforward to analyze, even in the delicate two- 
limit case n oo, m — associated with a metal. As mentioned at the 
beginning of Sec. 3.1 that it would seem most natural to analyze the limit 
of a metal assuming the standard Casimir configuration of parallel plates. 
In some sense the situation seems in fact to be the reverse: the spherical 
geometry is more easy to analyze in the metallic limit than the planar one. 
The reason for this is obvious: once plane plates are involved, one becomes 
confronted with two infinite spatial dimensions (the linear extensions of the 
plates), which lead to mathematically more delicate issues. Recent investi- 
gations of the Casimir effect for perfectly conducting plates have been given 
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by Klimchitskaya, Mostepanenko, and Geyer p6| . 

2. The basic expression for the free energy F, Eq. (1), holds for arbitrary 
temperatures as well as for arbitrary (frequency) dispersion relations. In the 
special case of a real metal we find, when adopting the Drude relation (in our 
opinion the preferable one at low frequencies) that the m = TE mode does 
not contribute. The total m = free energy for a metal becomes accordingly 
one half of the conventional expression (8). 

3. For a dielectric (finite n) there is no limiting problem: the m = case 
does not contribute to the TE mode at all. Figure 7 shows the magnitude of 
the m = free energy (thus associated with the TM mode) relative to the 
full free energy F. The relative contribution from the m = term is seen 
to increase with temperature, as one would expect physically; the relative 
weight of the low Matsubara frequencies becoming enhanced at high T. 

4. It ought to be stressed that F means everywhere the mutual free 
energy between the two concentric dielectric bodies. Thus F — when 
d = [b — a) ^ oo. It may be of interest to calculate the mutual internal 
energy E also. By means of the thermodynamic relation E = d{(3F)/d(3 we 
find immediately from Eq. (1) 



^ = -E E(2^ + i) 

m=0 1=1 



dxp' ^ 1 dxr 



XP' d(3 l-Xr d/3 



(27) 



Here the partial derivatives with respect to (3 are most conveniently calcu- 
lated on an analytic computer, on the basis of the expressions (2)- (6). If series 
approximation for the eigenvalues Xf^' or Xf^ were accessible, for instance 
in either of the temperature limits, it would be convenient to use Eq. (27) for 
evaluating approximate expressions for E. Obviously, E —>■ when d oo. 

5. It should be noted that since the spherical two-surface geometry that 
we are considering in our paper is different from the conventional parallel- 
plates geometry, this becomes reflected in the way in which we define the 
nondimensional temperature: we define it as t = 2iTa/f3, i.e., with the inner 
radius a as the geometrical variable instead of the conventional gap distance 
d. This is a natural definition in the case of curved surfaces. There are thus 
two different temperature scales involved here. This implies that at ordi- 
nary room temperature our problem becomes a high-temperature problem: 
by taking T = 300 K, a = 1 mm we obtain t to be as large as about 830. 
Under these circumstances, it follows from our Figs. 6 and 7 that the m = 
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TE term would be most important. A measurement of the surface surface 
force in this case would thus be critical. So far, no measurement of this force 
exists, however. So far, to our knowledge no conflict between our theory and 
experiments has been found. 

6. Recently, it has been argued by Klimchitskaya and Mostepanenko [|17 



and Bordag et al. |]T6| that the Drude dispersive model leads to inconsisten- 



cies at low frequencies, even in the conventional case of planar geometry. The 
reason for this, according to these authors, is that the Drude relation leads 
to a discontinuity in the reflection coefficient r2 as the imaginary frequency 
lu) — i> 0, in the case of perpendicular polarization. The plasma dispersion 
relation, instead of the Drude relation, is accordingly given preference by 
these authors since this discontinuity is not found to be present if one uses 
the plasma relation. 

These arguments are quite interesting, since they raise doubts not only 
about the validity of the Drude relation as such, but more generally even 
about the applicability of the Lifshitz formula at low temperatures. We 
intend to return to a study of this problem in a later paper . The problem 
with a) — * is most naturally discussed if one assumes planar geometry from 
the outset. Here, let us merely make a few remarks, related to our treatment 
above, choosing a specific metal for concreteness. When the radii a and h 
are large, our spherical system of course approaches that of planar geometry. 
Let us take aluminum, for which one has [|17|, ^ 

ujp = 1.9 X 10^S-\ 

7 = 9.6 X lO^^s"^ (28) 

For parallel plates separated by a gap d, the Matsubara frequencies (in 
dimensional units) are Um = 2TTkBTm/h. Zero-temperature theory is appli- 
cable as long as dkBT/hc ^1. At T = 0, the contribution from m = 
is negligible, since this contribution is completely buried in the Matsubara 
frequency integral. 

Assume now room temperature, T = 300 K. Then for aluminum 

cl;™ = (2.48m) X 10^S-\ (29) 

which shows that the difference between two adjacent Matsubara frequen- 
cies is in this case quite appreciable. The most important frequencies for 
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the Casimir effect occur when ujm ~ I'ncjd^ corresponding to the ordinary 
frequency uJm/2TT being of the order of the inverse transit time for photons 
between the two surfaces. This corresponds to m ~ hc/{dkBT). Taking for 
definiteness the gap tohe d = 0.5 fim, we obtain m ~ 15 to be the most sig- 
nificant Matsubara numbers. This is so far separated from m = that one 
should without any further calculation expect the contribution from m = 
to be quite small. And this agrees with the about 1 % level of temperature 



correction following from a more detailed calculation [|T^ . 

It is instructive to calculate also the conductivities, and the refractive 
indices, as following from the Drude model for the two lowest frequencies. For 
convenience we now use SI units. We first write the square of the refractive 
index, = e/eo, in the same form as in conventional low- frequency theory 
for metals: 

n\^u;) = l + ^. (30) 

EqUJ 

Here a{iuj) is an effective frequency-dependent conductivity. The Drude 
model, Eq. (24), corresponds to 

a{^u:) = (31) 

For the static conductivity, using Eq. (1), we find 

a(0) = 3.33 X lO^S/m, (32) 

whereas for the m = 1 case 

(T(icl;i) = 0.93 X lO^S/m. (33) 

The effective conductivity thus diminishes quickly when we move away from 
m = 0. The corresponding square refractive indices are 

n^{iu ^ 0) = (3.76/^:^) x 10^^ (34) 



n^{iui) = (1.05/cJi) X 10^^ (35) 

These refractive indices are large. However, the important point is that 
when a; — »• 0, Eq. (34) shows explicitly how n{iu})i2} when a; 0. Thus, 
this approach is in agreement with our option B above, in Sec. 3. 
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Let us return to the reflection coefficient r2 mentioned above, for perpen- 
dicular polarization. It is defined as [|17| 



r2 = (36) 

p + s 

where s and p are the conventional Lifshitz variables for planar geometry: 

k^^\k^\ = {Co/c){p'-iy/'. (37) 

The important question is: does r2 really becomes discontinuous at ci) = 
if one uses the Drude model? In our opinion, it does not. This can be seen 
from a power expansion in of the expressions above (we keep k± fixed; 
any normal metal must have a finite relaxation frequency 7). To lowest order 
we obtain s — p ^ u!p/{2'jk±c), s + p ^ 2k_ic/u, resulting in 

- 7^-- (38) 

This shows that r2 goes to zero smoothly (in our case linearly) as c!)/7 — ^ 0; 
no singularity at a; = is found. 

We intend to discuss these points in more detail in the mentioned forth- 
coming paper [^]. There, we will also discuss the recent claim of Fischbach 
et al. [29 1 that the results of Bostrom and Sernelius come into conflict 
with experiment as well as with basic thermodynamics. 

7. Generally, when comparing the outcome of Casimir calculations with 
experiments, care should be taken if the calculation involves summation over 
infinite series. It should here be observed that our discussion on the m = 
TE term in Section 3, as well as in our previous paper [Q, was based upon 
statistical methods, thus not involving summation methods like the one used 



by van Kampen et al. and others. Our viewpoint is quite physical: the 



static mode has to occur only once, not twice; it corresponds to the electric 
field being directed radially, thus transversely to the two spherical surfaces. 
This is precisely the TM mode. 

8. A remark on the so-called proximity force hypothesis [^] is in order, as 



this hypothesis is being made use of in connection with Casimir calculations 
for test bodies having spherical segments. Some doubts have been expressed 
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in the literature concerning the accuracy of this hypothesis. The issue has 
recently been analyzed by Barton (personal communication), with the result 
that this hypothesis remains valid in fourth-order as well as in second-order 
perturbation theory. Actually the hypothesis holds to all orders, as was 
shown by Langbein 



9. Finally, the recent experiment of Bressi et al. is interesting, since 
it reports a measurement of the Casimir force between conducting surfaces in 
a parallel configuration. At present, an accuracy at a 15 % level is achieved. 
It is to be hoped that this accuracy can be improved, although a direct 
experiment of this sort is obviously quite demanding. Again, as this is a 
low-temperature experiment, we cannot at present decide whether it is in 
accordance with our theory or not. 
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Figure Captions 



Figure 1 Sketch of the spherical geometry. 

Figure 2 Logarithm of mutual nondimensional free energy, \og^Q{— (3 Ft), 
versus relative width d/a for various values of the nondimensional tempera- 
ture t = 2'Ka/ (3. Refractive index n = 1.1. 

Figure 3 Logarithm log-^^Q{—f3F) versus d/a, when n = 1.1. 

Figure 4 Variation of logiQ^—PFt) versus logi^t for various values of 
d/a, when n = 1.1. 

Figure 5 Same as Fig. 4, but with n = 2.0. 

Figure 6 Same as Fig. 5, but for an ideal metal (n = cx), the m = 
mode counted twice). Reproduced (with corrected labeling) from Ref. []1|. 

Figure 7 Relative importance of the zero frequency term m = in the 
free energy, versus t = 2TTa/ (3. The quantity Y is defined in Eq. (26). 
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